Search for: All records

Characterized by similar-to-today CO2 (∼400 ppm) and surface temperatures approximately 3°–4°C warmer than the preindustrial, the mid-Pliocene warm period (mPWP) has often been used as an analog for modern CO2-driven climate change and as a constraint on the equilibrium climate sensitivity (ECS). However, model intercomparison studies suggest that non-CO2boundary conditions—such as changes in ice sheets—contribute substantially to the higher global mean temperatures and strongly shape the pattern of sea surface warming during the mPWP. Here, we employ a set of CESM2 simulations to quantify mPWP effective radiative forcings, study the role of ocean circulation changes in shaping the patterns of sea surface temperatures, and calculate radiative feedbacks during the mPWP. We find that the non-CO2boundary conditions of the mPWP, enhanced by changes in ocean circulation, contributed to larger high-latitude warming and less-stabilizing feedbacks relative to those induced by CO2alone. Accounting for differences in feedbacks between the mPWP and the modern (greenhouse gas–driven) climate provides stronger constraints on the high-end of modern-day ECS. However, a quantification of the forcing of non-CO₂boundary condition changes combined with the distinct radiative feedbacks that they induce suggests that Earth system sensitivity may be higher than previously estimated. 

This dataset originates from a new CESM2 CAM6 perturbed parameter ensemble (PPE) designed to explore climate and hydroclimate dynamics under a wide range of sea surface temperature (SST) conditions. The SST varies from 4 degrees Celsius colder to 16 degrees Celsius warmer than preindustrial levels, encompassing a broad spectrum of mean temperatures spanning the past 65 million years. This dataset offers valuable insights into climate and hydroclimate responses, as well as weather and climate extremes under diverse conditions.The dataset includes results from nine PPE simulations with different SST scenarios: preindustrial (PREI), 4K cooler (M04K), and 4K, 8K, 12K, and 16K warmer (P04K to P16K). For SSTs exceeding 8K warming, sea ice was removed to improve numerical stability. Each PPE set consists of 250 ensemble members, with 45 parameters related to microphysics, convection, turbulence, and aerosols perturbed using Latin Hypercube Sampling. An additional simulation with default parameter settings brings the total to 251 simulations, each running for five years using CAM6.3 (https://github.com/ESCOMP/CAM/tree/cam6_3_026; with additional paleo modifications).Post-processing converted the data into compressed NetCDF-4 format. All 251 runs were concatenated using ncecat to minimize the number of files. For example, the following file contains monthly surface temperature data from the preindustrial PPE: f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251/atm/proc/tseries/month_1/f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251.cam.h0.TS.000101-000512.ncA detailed variable list [https://rda.ucar.edu/OS/web/datasets/d651038/docs/detailed_vars.txt] can be found in the Documentation Tab.Parameter values are provided in the PPE Parameter File. More details can be found in the paper: Zhu et al. (2025). Investigating the State Dependence of Cloud Feedback Using a Suite of Perturbed Parameter Ensembles, Journal of Climate. 

Past climate states hold valuable insights into future climate change. Among those states, mid-Pliocene (3.3 - 3.0 Ma) is often studied as an important analog to near future climate change following an intermediate warming pathway. This time interval featured topography and geography like present-day, yet with retreated polar ice sheets and expanded boreal forests, potentially reflecting equilibrium earth system responses to CO2 forcing at a centennial to millennial time scale. Despite the prolific research on Pliocene climate, little is known about the amount of radiative forcing, especially from changing boundary conditions, that drives the Pliocene climate. Existing constraints mainly focused on well-mixed greenhouse gases and aerosols. Here, we applied the methodology commonly used to quantify radiative forcing of future climate and its sources to constrain radiative forcing of the mid-Pliocene climate using three generations of Community Earth System Models (CCSM4, CESM1.2, and CESM2). To calculate ERF, the difference in net top of the atmosphere radiative fluxes is computed between a pre-industrial control and a mid-Pliocene simulation. Both are carried out with prescribed pre-industrial sea surface temperature. The three mid-Pliocene simulations separately feature a 400 ppm CO2 (the level of mid-Pliocene), mid-Pliocene geography and topography, and mid-Pliocene ice and vegetation. Changing atmospheric temperature, water vapor, surface albedo, and clear vs total sky radiative fluxes are further extracted from these simulations to calculate radiative adjustments with published radiative kernels for CESM. In our preliminary results with CESM1.2, we found that ERF is 1.754 W m-2 for CO2 forcings, 1.143 W m-2 for vegetation and ice sheet forcing, and -0.339 W m-2 for geographic and topographic forcing. Further, ERF from boundary condition changes mostly arises from changing surface albedo with 1.626 W m-2 for vegetation and ice sheet changes and –0.54 W m-2 for geographic and topographic changes respectively. Radiative adjustments from water vapor responses tend to amplify the instantaneous forcing with the most profound effect induced by vegetation and ice sheet changes. These results underscore the importance of constraining radiative forcing from changes in boundary conditions, which is potentially key to understanding drivers of past climate warmth and inter-model spread in simulated past climate states. 

The cross-equatorial southwesterly winds from the eastern equatorial Pacific direct moisture toward the Pacific coast of northwestern South America, where subsequent orographic lifting creates the wettest regions in the world. The Choco low-level jet is emblematic of broader westerly winds in this region and is projected to weaken by the end of the 21st century, but climate models show considerable disagreement about the extent of weakening. Using contemporary observations, we demonstrate that the configuration of westerly winds in the eastern equatorial Pacific is reflected by hydrogen isotopes in precipitation (δDp) in western Ecuador. As westerly winds strengthen, δDp increases from greater transport of δDvapor enriched in deuterium from the Eastern Pacific Warm Pool. We apply this framework to a new record of reconstructed δDp using leaf waxes in ocean sediments off the coast of Ecuador (ODP1239, 0◦40.32′ S, 82◦4.86′ W) that span the Plio-Pleistocene. Low δDp in the early Pliocene indicates weak westerly water vapor transport in a warmer climate state, which is attributed to a low sea surface temperature gradient between the cold tongue and off-equatorial regions in the eastern equatorial Pacific. Near 3 Ma, westerly water vapor transport weakens, possibly as a result of shifts in the Intertropical Convergence Zone forced by high latitude Northern Hemisphere cooling. In complementary isotope-enabled climate simulations, a weak Choco jet and westerly water vapor transport in the early Pliocene are matched by a decrease in δDp and hydroclimate changes in western Ecuador. Precipitation from the Choco jet can cause deadly landslides and weakened westerly winds in the early Pliocene implies a southward shift of these hazards along the Pacific coast of northwestern South America in the future. 

Abstract The state dependence of cloud feedback—its variation with the mean state climate—has been found in many paleoclimate and contemporary climate simulations. Previous results have shown inconsistencies in the sign, magnitude, and underlying mechanisms of state dependence. To address this, we utilize a perturbed parameter ensemble (PPE) approach with fixed sea surface temperature (SST) in the Community Atmosphere Model, version 6. Our suites of PPEs span a wide range of global mean surface temperatures (GMSTs), with spatially uniform SST perturbations of −4, 0, 4, 8, 12, and 16 K from the preindustrial. The results reveal a nonmonotonic variation with GMSTs: Cloud feedback increases under both cooler and warmer-than-preindustrial conditions, with a rise of ∼0.1 W m⁻²K⁻¹under a 4-K colder climate and ∼0.4 W m⁻²K⁻¹under a 12-K warmer climate. This complexity arises from differing cloud feedback responses in high and low latitudes. In high latitudes, cloud feedback consistently rises with warming, likely driven by a moist adiabatic mechanism that influences cloud liquid water. The low-latitude feedback increases under both cooler and warmer conditions, likely influenced by changes in the lower-tropospheric stability. This stability shift is tied to nonlinearity in thermodynamic responses, particularly in the tropical latent heating, alongside potential state-dependent changes in tropical circulations. Under warmer-than-preindustrial conditions, the increase in cloud feedback with warming is negatively correlated with its preindustrial value. Our PPE approach takes the model parameter uncertainty into account and emphasizes the critical role of state dependence in understanding past and predicting future climates. Significance StatementThis study focuses on how cloud feedback—one of the most uncertain aspects of climate change—varies as global temperatures rise. We found that the cloud feedback decreases at first with warming and then increases, showing significant variation. This complexity stems from nonlinear thermodynamics, such as the Clapeyron–Clausius relationship, which describes how temperature affects moisture in the atmosphere. Our results indicate that the cloud feedback depends on the level of global warming, which is a significant factor rooted in fundamental physics. Recognizing this dependence is important for studies that aim to interpret past climates and predict future climate changes. 

Abstract. The mid-Pliocene Warm Period (mPWP, 3.3–3.0 Ma) was characterised by an atmospheric CO2 concentration exceeding 400 ppmv with minor changes in continental and orbital configurations. Simulations of this past climate state have improved with newer models but still show some substantial differences from proxy reconstructions. There is little information about atmospheric aerosol concentrations during the Pliocene, but previous work suggests that it could have been quite different from the modern period. Here we apply idealised aerosol scenario experiments to examine the importance of aerosol forcing on mPWP tropical precipitation and the possibility of aerosol uncertainty explaining the mismatch between reconstructions and simulations. The absence of industrial pollutants leads to further warming, especially in the Northern Hemisphere. The Intertropical Convergence Zone (ITCZ) becomes narrower and stronger and shifts northward after removal of anthropogenic aerosols. Though not affecting the location of monsoon domain boundary, removal of anthropogenic aerosol alters the amount of rainfall within the domain, increasing summer rain rate over eastern and southern Asia and western Africa. This work demonstrates that uncertainty in aerosol forcing could be the dominant driver in tropical precipitation changes during the mid-Pliocene: causing larger impacts than the changes in topography and greenhouse gases. 

Equilibrium climate sensitivity (ECS) quantifies the amount of warming resulting from a doubling of the atmospheric CO2 forcing. Despite recent advancements in climate simulation capabilities and global observations, there remains large uncertainty on the degree of future warming. To help alleviate this uncertainty, past climates provide a valuable insight into how the Earth will respond to elevated atmospheric CO2. However, there is evidence to suggest that ECS is dependent on background climate warmth, which may interfere with the direct utilization of paleo-ECS to understand present-day ECS. Thus, it is important that a range of different climate states are considered to better understand the factors modulating the relationship between CO2 and temperature. In this study, we focus on three time intervals: the mid-Pliocene Warm Period (3.3 – 3.0 Ma), the mid-Miocene (16.75 – 14.5 Ma), and the early Eocene (~50 Ma), in order to sample ECS from Cenozoic coolhouse to hothouse climates. Here, we combine the Bayesian framework of constraining the ECS and its uncertainty with several published methods to estimate the global mean surface temperature (GMST) from sparse proxy records. This framework utilizes an emergent constraint between the simulated GMST changes and climate sensitivities across the model ensemble. For each time interval, we employ a combination of parametric and non-parametric functions, coupled with a probabilistic approach to derive a refined estimate. Preliminary results for the Pliocene indicate a GMST reconstruction of approximately 19.3°C, which is higher than previous estimates that were derived using only marine records. Using this estimate, we calculate an ECS that is also higher than previously published values, especially due to the inclusion of high-latitude terrestrial temperature records into our estimates. Intriguingly, using the consistent methodology, our calculated ECS for the early Eocene is lower than that of the mid-Pliocene. This result does not support an amplified ECS in hothouse climate, and points to a potentially important role of ice albedo feedback in amplifying the ECS in coolhouse climate. Ongoing work will apply the same methodology to the mid-Miocene and further investigate the source for the estimated ECS state dependency between these climate intervals. 

Southwestern North America is currently experiencing a multidecadal megadrought, with severe consequences for water resources. However, significant uncertainty remains about 21st century precipitation changes in this semi-arid region. Paleoclimatic records are essential for both contextualizing current change, and for helping constrain the sensitivity of regional hydroclimate to large-scale global climate. In this paper, we present a new 2.8 Ma to present compound-specific isotopic record from Clayton Valley, the site of a long-lived paleolake in the southern Great Basin. Hydrogen and carbon isotopes from terrestrial plant leaf waxes provide evidence of past shifts in rainfall seasonality as well as ecosystem structure, and help contextualize the formation of this lithium-rich lacustrine basin. Our results suggest that regional hydroclimates underwent a substantial reorganization at the Plio-Pleistocene boundary, especially between 2.6 and 2.0 Ma. In this interval, a reduced latitudinal temperature gradient in the North Pacific likely resulted in a northward shift in storm tracks, and a reduction in winter rainfall over the southern Great Basin. This occurred against a background of increased summer rainfall and a greater accumulation of lithium in the lake basin. Our interpretation is corroborated by a compilation of Plio-Pleistocene north Pacific sea surface temperature records, as well as an isotope-enabled model simulation. Overall, these results suggest that past shifts in rainfall seasonality helped set the stage for the development and dessication of lithium-rich lacustrine deposits. 

Abstract In August 2022, Death Valley, the driest place in North America, experienced record flooding from summertime rainfall associated with the North American monsoon (NAM). Given the socioeconomic cost of these type of events, there is a dire need to understand their drivers and future statistics. Existing theory predicts that increases in the intensity of precipitation is a robust response to anthropogenic warming. Paleoclimatic evidence suggests that northeast Pacific (NEP) sea surface temperature (SST) variability could further intensify summertime NAM rainfall over the desert southwest. Drawing on this paleoclimatic evidence, we use historical observations and reanalyzes to test the hypothesis that warm SSTs on the southern California margin are linked to more frequent extreme precipitation events in the NAM domain. We find that summers with above-average coastal SSTs are more favorable to moist convection in the northern edge of the NAM domain (southern California, Arizona, New Mexico, and the southern Great Basin). This is because warmer SSTs drive circulation changes that increase moisture flux into the desert southwest, driving more frequent precipitation extremes and increases in seasonal rainfall totals. These results, which are robust across observational products, establish a linkage between marine and terrestrial extremes, since summers with anomalously warm SSTs on the California margin have been linked to seasonal or multi-year NEP marine heatwaves. However, current generation earth system models (ESMs) struggle to reproduce the observed relationship between coastal SSTs and NAM precipitation. Across models, there is a strong negative relationship between the magnitude of an ESM’s warm SST bias on the California margin and its skill at reproducing the correlation with desert southwest rainfall. Given persistent NEP SST biases in ESMs, our results suggest that efforts to improve representation of climatological SSTs are crucial for accurately predicting future changes in hydroclimate extremes in the desert southwest.